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PERSP I RATION  POISONING  OF  PROTECTIVE  CLOTHING  MATERIALS 

PART  II 

MATHEMATICAL  MODEL  FOR  A  COMPLEX  ADSORPTION  BED 
INTRODUCTION 

The  mathematical  modeling  of  adsorption  processes  is  well  documented 
and  the  literature  on  the  subject  is  extensive.  One  of  the  most  complete 
and  useful  approaches  is  that  proposed  by  Schneider  and  Smith^  for  a  bed 
of  particles  in  which  adsorption  takes  place  as  a  three-step  mechanism: 

(1)  diffusion  from  the  bulk  gas  phase  to  the  external  surface  of  the  par¬ 
ticle,  (2)  diffusion  into  the  particle,  and  (3)  adsorption  on  the  particle's 

2 

surface.  Masamune  and  Smith  have  presented  models  that  describe  the  above 
mechanisms  as  controlling  resistances,  both  singly  and  in  pairs,  and  offer 
an  alternate  approach  to  modeling  an  adsorption  process.  If  an  adsorption 
process  can  be  characterized  by  a  single  controlling  resistance,  the  mathe¬ 
matical  work  and  subsequent  model  are  greatly  simplified.  This  in  turn  gives 
the  investigator  a  better  understanding  of  the  ways  in  which  the  physical 
properties  of  his  adsorption  system  affect  the  adsorption  process.  The  util¬ 
ity  of  such  models  for  well-defined  systems  has  been  demonstrated  by  the  ex¬ 
cellent  agreement  of  experimental  and  theoretical  transmission  curves  of 

.  1  2 

C/CQ  vs.  time  obtained  by  Schneider  and  Smith  and  Masamune  and  Smith.  How¬ 
ever,  problems  arise  in  this  procedure  when  the  adsorption  system  is  complex. 


Schneider,  P.,  and  J.  M.  Smith,  AIChE  J. ,  14,  762  0968). 
^Masamune,  S.,  and  J.  M.  Smith,  AIChE  J. ,  11,  34  (1965). 


not  well-defined,  and  of  questionable  homogeneity.  The  work  performed  dur¬ 
ing  this  study  was  on  such  a  system  in  which  the  adsorbent  was  activated 
carbon  impregnated  in  foam  material  with  a  nylon  backing.  Inherent  uncer¬ 
tainties  for  this  system  included  the  amount  of  carbon  in  a  foam  sample, 
carbon  particle  size,  non- homogeneity  of  foam  material  samples,  and  the 
characterization  of  flow  through  a  foam  matrix  embedded  with  charcoal. 

The  adsorbent  material  was  developed  by  the  U.  S.  Army  Natick  R&D 
Command  for  use  as  a  protective  overgarment  to  adsorb  chemical  agents, 
especially  poisonous  gases.  The  material  was  found  very  effective  for 
short  periods  of  time,  but  due  to  its  thermal  insulating  properties  its 
adsorptive  capacity  for  toxic  gases  was  reduced  because  of  "poisoning" 
human  perspiration.  The  long-time  objective  of  the  grant  that  funded  this 
research  is  to  remedy  this  "poisoning"  problem;  therefore,  the  modeling 
work  considersthe  different  types  of  conditions  described  in  the  experi¬ 
mental  procedure. 

The  material  consisted  of  a  layer  of  polyurethane  foam  bonded  to  a  nylon 
tricot  and  impregnated  with  activated  carbon  held  on  the  material  with  a 
polymer  latex  binder.  Thickness  of  the  material  was  approximately  0.18  cm. 
Samples  of  this  carbon  impregnated  foam  were  received  as  bolts  of  material 
which  were  found  to  have  wide  variations  in  adsorptive  capacity.  Work  done 
during  this  study  was  on  the  second  and  third  bolts  received  from  Natick 
and  will  be  referred  to  as  bolt  2  and  bolt  3,  respectively.  A  detailed 
description  of  the  material  is  given  In  Part  I  of  this  report. 

The  Initial  model  work  was  directed  toward  developing  the  three  para¬ 
meter  model  of  Schneider  and  Smith^  and  applying  it  to  the  adsorption  of 


op.  cit. 
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carbon  tetrachloride  vapor  by  charcoal  impregnated  foam.  It  was  decided 
to  use  the  method  of  moments  to  evaluate  the  three  rate  parameters  from 
the  experimental  data.  This  procedure  required  three  independent  moment 
equations  and  consistent  data  with  a  minimum  of  error  magnification  for 
higher  moments.  The  zeroth  moment  equation  gave  an  expression  for  deter¬ 
mining  «A,  the  equilibrium  adsorption  constant,  and  the  first,  second, 
and  third  moment  equations  gave  independent  expressions  for  the  three  rate 
parameters.  The  solution  of  these  equations  yielded  negative  and  imaginary 
values  for  two  of  the  rate  parameters  and  revealed  the  inadequacy  of  the 
data  for  use  in  evaluating  more  than  the  first  moment  numerically.  This 
made  necessary  a  simplified  approach  of  assuming  a  single  resistance  rate 
controlling  and  evaluating  this  resistance  from  the  first  moment  equation. 
Diffusion  into  the  particle  was  assumed  to  be  the  rate  determining  mechanism 
and  the  model  equations  were  solved.  The  model  worked  well  in  predicting 
breakthrough  curves,  but  it  was  found  that  the  models  obtained  by  assuming 
either  external  diffusion  or  surface  adsorption  rate  controlling  also  fit 
the  data  well.  Differentiation  between  models  was  further  complicated  by 
an  inability  to  change  the  adsorbent's  physical  properties  and  by  a  limited 
range  of  temperature  and  flow  in  the  adsorption  apparatus. 

Subsequent  mathematical  work  showed  that  when  the  rate  parameter  was 
evaluated  from  the  first  moment,  the  external  diffusion  and  surface  adsorp¬ 
tion  models  were  exactly  the  same,  and  that  the  model  assuming  pore  diffu¬ 
sion  controlling,  while  not  mathematically  identical,  was  not  significantly 
different  from  the  other  two  cases.  It  was  then  proposed  that  the  first 
moment  equation  for  the  three  parameter  model  be  expressed  as  an  overall 
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coefficient,  UA: 

Hi  (Hl_  .  1)  *  r_L  *  I  +  _J i  =  ± 

Q  'M  2  2'  ‘■kfA  5DCA/R  kads6J  UA 

This  overall  coefficient  would  reflect  the  influence  of  the  three  different 
resistances  and  would  make  possible  the  use  of  a  model  with  only  one  rate 
parameter  when  the  controlling  mechanisms  of  mass  transfer  were  not  known. 

An  overall  coefficient  model  was  developed  and  when  the  rate  parameter  was 
evaluated  from  the  first  moment,  this  model  was  shown  to  be  the  same  as  the 
model  of  external  diffusion  rate  controlling  and  the  same  as  the  model  of 
surface  adsorption  rate  controlling.  It  was  also  found  that  the  overall 
coefficient  model  produced  results  that  were  not  significantly  different 
from  those  produced  by  the  pore  diffusion  controlling  model.  This  procedure 
offered  an  alternative  approach  in  modeling  an  adsorption  process  in  which 
the  adsorbent  properties  and  adsorption  conditions  could  not  be  varied  ef¬ 
fectively  to  determine  the  true  controlling  mechanisms  of  mass  transfer. 
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EXPERIMENTAL 

A  schematic  diagram  of  the  apparatus  for  the  study  of  the  dynamics 
of  carbon  tetrachloride  vapor  adsorption  by  charcoal  impregnated  foam 
material  is  shown  in  Figure  1  with  a  list  of  components  in  Table  1.  This 
apparatus  consisted  of  a  flow  system  of  stainless  steel  tubing  with  an 
adsorption  chamber  enclosed  in  a  plexiglass  box  for  constant  temperature 
control.  A  pure  nitrogen  stream  from  the  nitrogen  supply  cylinder  was 
split  at  the  inlet  and  sent  through  lines  (1)  and  (2).  Nitrogen  in  line 
(2)  flowed  through  the  carbon  tetrachloride  bubbler  (4)  immersed  in  an 
ice-bath  to  saturate  the  nitrogen  stream  with  carbon  tetrachloride  vapor 
at  0°C.  The  pure  nitrogen  stream  in  line  (1)  flowed  through  an  orifice 
(5)  with  a  hook-gage  manometer  (6)  and  rotameter  in  series  and  was  mixed 
with  the  nitrogen-carbon  tetrachloride  stream  for  dilution  to  the  desired 
concentration.  The  flow  control  of  the  pure  nitrogen  stream  in  line  (1) 
was  critical  because  of  the  high  dilution  factor.  This  flow  rate  wa.,  set 
roughly  by  the  rotameter  (8)  and  then  adjusted  by  measuring  the  pressure 
drop  across  the  orifice  to  one  ten-thousandth  of  an  inch  of  water  with  the 
hook-gage  manometer. 

The  diluted  nitrogen-carbon  tetrachloride  mixture  was  sent  into  the 
constant  temperature  box  (25)  through  a  temperature  equilibrating  coil  (14) 
and  into  a  manifold  (17).  A  vapor  stream  was  drawn  from  the  manifold  (at 
18)  and  sent  through  the  sample  holder  (19)  at  a  flow  rate  set  by  rotameter 
(22).  This  sample  holder  consisted  of  two  stainless  steel  cups  between 
which  the  foam  material  was  sandwiched,  tightened  with  a  clamp,  and  sealed 
with  wax.  The  top  cup  contained  a  perforated  metal  sheet  to  assure  uniform 
gas  flow  through  the  sample.  The  vapor  stream  leaving  the  sample  holder 
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was  sent  through  a  gas  collection  coil  on  the  chromatograph  and  the  remainder 
of  the  vapor  mixture  entering  the  manifold  was  vented  through  an  exhaust 
hood  (23). 

Analysis  of  the  carbon  tetrachloride  concentration  in  the  stream  exit¬ 
ing  the  sample  cup  was  made  by  a  Perkin-Elmer  gas  chromatograph  using  a 
column  of  silicone  oil  D.C.  #200  with  a  thermal  conductivity  detector.  A 
Moseley  strip  chart  recorder  monitored  the  chromatograph  output.  Reference 
peaks  for  the  initial  concentration  were  obtained  at  the  beginning  and  end  of 
each  adsorption  run  by  taking  a  sample  stream  from  the  manifold  at  (26)  and 
sending  it  through  the  chromatograph.  Injections  were  made  using  the  gas 
sampling  valve  of  the  chromatograph  at  three  minute  intervals,  and  from  the 
recorded  output  (chromatographic  curves),  transmission  curves  of  C(t)/CQ  vs. 
time  were  generated. 

A  method  of  conditioning  the  foam  samples  was  devised  to  improve  con¬ 
sistency  from  sample  to  sample.  The  foam  material  was  cut  into  five  inch 
circles,  soaked  in  distilled  water  for  twenty-four  hours,  wrung  out  between 
rubber  rollers,  and  then  allowed  to  equilibrate  in  a  room  maintained  at  70°F 
(21°C)  and  65%  relative  humidity.  Treatment  with  other  solutions  (e.g., 
sweat,  lactic  acid)  was  accomplished  by  following  the  water  conditioning 
with  treatment  with  the  proper  solution  and  again  allowing  the  sample  to 
equilibrate  at  70°F  (21°C)  and  65%  relative  humidity.  Thu^  valid  compari¬ 
sons  could  then  be  made  between  different  treatments  and  also  among  condi¬ 
tioned  samples  for  different  run  conditions. 

The  adsorption  bed  was  composed  of  either  single  or  multiple  layers 
of  the  carbon  impregnated  foam  material:  for  this  study  beds  of  one  and 
three  layers  were  considered.  In  this  manner  the  effects  of  different  bed 
depths  could  be  studied. 
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Carbon  tetrachloride  was  suggested  as  the  adsorbable  vapor  since  cor¬ 
relations  exist  between  carbon  tetrachloride  and  the  toxic  gases  in  which 
the  Army  was  interested. 
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DISCUSSION 


One  of  the  most  useful  and  complete  approaches  to  modeling  an  adsorp¬ 
tion  process  is  that  of  Schneider  and  Smith^  who  proposed  that  the  adsorp¬ 
tion  of  a  gas  flowing  through  a  bed  of  spherical  particles  could  be  modeled 
as  a  three  step  mechanism:  diffusion  from  the  bulk  gas  to  the  external  sur¬ 
face  of  the  particle,  pore  diffusion  into  the  particle,  and  adsorption  on 
the  particles  surface.  For  the  present  study  the  adsorption  bed  was  com¬ 
posed  of  carbon  particles,  assumed  to  be  spherical,  embedded  in  an  inert 
matrix  of  urethane  foam.  Since  very  thin  beds  were  used,  the  effect  of  axial 
dispersion  was  neglected  and  plug  flow  was  assumed.  The  presence  of  inerts 
in  the  adsorption  bed  necessitated  defining  two  void  fractions,  a  and  e, 
where  a  represents  void  volume/total  volume  that  must  be  used  in  determining 
the  volume  in  which  convection  and  accumulation  take  place  in  the  gas  phase, 
and  where  1  -  e  represents  carbon  volume/total  volume  which  must  be  used  in 
expressing  the  rate  of  removal  in  terms  of  the  volume  of  carbon  particles. 
With  these  modifications  a  mass  balance  of  the  adsorbable  component  in  the 
gas  phase  yields: 
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The  rate  of  adsorption  (assumed  to  be  linear)  is: 
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with  the  various  boundary  conditions  listed  below: 
external  diffusion  boundary  condition: 


DC  (1F>r=R  =  kflC  "  ^ 


internal  diffusion  boundary  condition: 


=  0  at  r  =  0  for  t  >  0 


initial  conditions: 

C  =  0  at  z  >  0  for  t  =  0  (6) 

C.  =  0  at  r  >  0  for  t  =  0  (7) 

C  =  CQ  at  z  *  0  for  t  >  0  (8) 

Equations  1  through  8  have  been  solved  by  Rosen“  using  Laplace  trans¬ 

forms  and  the  inversion  integral  to  give  an  expression  for  C(t),  the  gas 
concentration  exiting  the  bed,  in  the  form  of  an  infinite  integral.  To  use 
Rosen's  expression  for  C(t),  it  is  necessary  to  know  the  values  of  the  three 
rate  parameters,  k^,  kacjs ,  and  Dc,  for  the  adsorption  system  of  this  study. 

An  empirical  approach  to  evaluating  these  parameters  was  ruled  out  due  to  un¬ 
certainties  in  the  adsorbent's  properties  and  shortcomings  in  the  adsorption 
apparatus.  These  problems  include  (1)  an  inability  to  change  the  adsorbent's 
physical  properties,  such  as  particle  size  and  distribution,  (2)  an  uncertainty 
as  to  the  characterization  of  flow  through  a  foam  matrix  embedded  with  carbon, 
and  (3)  a  limited  range  of  flow  rate  and  temperature  in  the  adsorption  apparatus. 

The  method  of  moments  offered  an  alternative  approach  in  that  the  three 
rate  parameters  could  be  determined  analytically  for  each  adsorption  run  by 


3Rosen,  J.  B. .  J. 


Es_,  20,  387  (1952). 
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evaluating  the  first,  second  and  third  moments  numerically  and  using  the 

corresponding  moment  equations  to  solve  directly  for  k^,  kfids,  and  Dc. 

4 

Moment  equations  were  derived  for  a  step  input  using  Aris1  theorem  and  were 
found  to  be: 
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18 


Q°  (kfA)3  175(DCA/R)3 


(12) 


144 


Ckads6)3  5(kfA)2(0cA/R) 


+  +  2§5)  ( 


27 


(kfA)(DcA/R) 


1 


18 


36 

Z 


18 


(kfA)  (kads6)  35(0  A/R)  (kads6)  (k  A)(kads6) 


5(DCA/R)(kad$6)2  5(DcA/R)(kfA)(kads6) 

For  this  system  in  which  constants  such  as  p^,  a,  e,  and  R  are  not  known 
and  are  very  difficult  to  measure  accurately,  constants  were  grouped  together 
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and  lumped  with  the  rate  parameter  to  give  a  coefficient  that  would  be 
representative  of  the  rate  parameter  contained  within  the  grouping.  An 
additional  advantage  in  this  procedure  is  the  conversion  from  a  weight 
basis  of  carbon  particles  (which  is  not  known)  to  a  weight  basis  of  foam 
material  by  assuming  that  6,  grams  carbon/gram  foam,  is  constant  for  any 
given  bolt  of  carbon  impregnated  foam  material.  The  following  substitu¬ 
tions  were  made: 

3 

If  =  grams  foam  material/cm  total  volume 
3  3 

1  -  e  =  cm  carbon/cm  total  volume 
3  3 

a  =  cm  void/ cm  total  volume 

2 

Af  =  bed  cross  section,  cm 

z  =  bed  depth,  cm 

6  =  g  carbon/g  foam 

3 

Ka'  =  equilibrium  adsorption  constant,  cm  adsorbed/g  carbon 


Then  the  weight  of  the  sample  bed,  W^,  is: 


Wt  =  AfZpm 


(13) 


and  the  volumetric  flow  rate  through  the  sample,  Q,  is: 


Q  =  AfVa  (14) 

So  that  Wt  z  pm  (15) 

Q  "  V  a 

3 

A  new  equilibrium  constant  K^,  cm  adsorbed/gram  foam,  can  be  defined  on  a 
weight  basis  of  foam  material  as: 

KA  =  6V  (16) 


The  apparent  foam  density  Is  shown  to  be: 


_  /,  \  cm  carbon _  1  g  foam  g  carbon 

%  cn,3  total  volume  '  s  9  carl,on  '  “p  cm3  carbon 


(17) 


1  - 


q  foam 


6  p  cm^  total  volume 


Substituting  (15),  (16),  and  (17)  into  the  zeroth  moment  equation: 

W. 


Mo  =  7 


Z  1  -  e  „  ,  "t 


Pp  V  =  t  KA 


08) 


The  quantity  A,  used  in  the  lumped  parameters  of  the  first,  second,  and 
third  moment  equations,  is  defined  as  surface  area  of  carbon  particles/unit 
weight  of  foam  material  and  is  shown  to  be: 


2  2  3 

^  _  4nR  cm  surface  area  .  1  cm  carbon  .  „  g  carbon  =  14  09) 


4/3  irR  cm  carbon 


p.~g  carbon 


g  foam 


The  substitution  of  (18)  and  (19)  into  the  higher  moment  equations  gave  the 
results  listed  as  equations  (10),  (11),  and  (12). 

Equations  for  the  first,  second,  and  third  moments  allowed  the  three 
rate  parameters  to  be  calculated  directly  from  a  single  adsorption  run;  how¬ 
ever,  problems  arose  because  of  the  generation  of  imaginary  or  negative 
values  for  k^A  and  kaj$.  This  was  probably  due  to  inaccuracies  in  data 
since  the  calculation  of  the  second  and  higher  moments  necessitates  very 
accurate  experimental  data.  Serious  troubles  arise  from  the  uncertainty 
of  the  exact  position  of  the  tails  of  the  transmission  curve  where  a  minor 
change  in  concentration  is  enlarged  by  the  second  (or  higher)  power  of  time. 
This  problem  led  to  the  development  of  a  model  in  which  one  of  the  three  re¬ 
sistances  is  considered  to  be  rate  determining  and  can  be  found  from  the 

2 

first  moment  equation.  Masamune  and  Smith  have  obtained  solutions  for  the 
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three  cases  where  each  resistance  is  considered  controlling  and  these  are 
listed  in  Table  2. 

The  kf  and  kflds  models  shown  in  Table  2  are  the  same,  and  when  the 
first  moment  equation  is  used  to  evaluate  the  controlling  resistance  as 
either  kfA,  kac|s,  or  DcA/R,  all  three  models  give  approximately  the  same 
results.  The  differentiation  between  models  was  further  complicated  by  an 
inability  to  change  the  adsorbent's  physical  properties  (such  as  particle 
size  and  distribution),  a  limited  range  of  inlet  concentration,  flow  and 
temperature  in  the  adsorption  apparatus,  and  as  previously  mentioned,  C/CQ 
transmission  curves  that  are  inadequate  for  complete  moments  analysis. 

An  alternative  way  of  viewing  this  adsorption  process  is  that  the  re¬ 
sistances  to  mass  transfer  from  the  bulk  gas  to  the  particle's  surface  may 
be  represented  by  one  overall  mass  transfer  coefficient.  The  form  of  the 
first  moment  equation  for  the  three  parameter  model  suggests  that  an  over¬ 
all  coefficient,  UA,  might  be  written  as: 


»'»t 


Mi  i 

(• \  - 1)  ■ 

V 


5D^7R 


The  three  resistances  to  adsorption  are  not  in  series  since  pore  dif¬ 
fusion  contributes  a  variable  resistance  that  is  a  function  of  the  spheri¬ 
cal  coordinate,  r,  and  may  occur  in  parallel  with  surface  adsorption.  How¬ 
ever,  one  may  picture  the  pore  diffusion  resistance  term,  1/5DCA/R,  as  an 
average  resistance  or  the  pore  diffusion  resistance  at  the  effective  path 
length  that  all  molecules  must  travel  before  being  adsorbed  on  the  surface. 
Then  the  definition  of  an  overall  coefficient  provides  a  means  of  accounting 
for  the  effects  of  the  three  resistances  as  an  average  or  overall  resistance. 
This  overall  coefficient  when  evaluated  from  the  first  moment  would  then 


i 


reflect  the  influences  of  each  of  the  three  resistances  and  would  simplify 
the  modeling  in  that  higher  moments  equations  would  be  unnecessary  in 
evaluating  UA.  The  development  of  an  overall  coefficient  model  is  out¬ 
lined  below: 


s  _  UA  / p  t/  i  p  \ 
at  "  WT  ltKA  "  Ladsj 


and  the  mass  balance  in  the  gas  phase  is: 

V  |£  +  +  Iz£  3Cads,  .  0 

az  at  Pp  a  at 

Derivation  of  the  moments  equations  for  the  above  model  yields: 

"o  ■  ^  ¥  »p  Ka'  *  V  KA 

(>VMo2  -  ,/2>  V1* '  m 

The  differential  equations  were  solved  using  the  procedure  outlined  by 

5 

Pigford  and  Marshall  and  were  also  put  into  the  form  presented  by  Masamune 
2 

and  Smith  for  comparison.  (See  Table  2.)  For  the  cases  of  kf  controlling 
or  kads  controlling  (i.e.,  making  the  substitutions  UA  =  kfA  or  UA  =  ka(js  6, 
respectively),  the  overall  coefficient  model  is  identical  to  the  two  correspond¬ 
ing  single  resistance  models.  For  the  case  of  pore  diffusion  controlling,  the 
overall  coefficient  model  gives  the  same  results  within  the  range  of  experi¬ 
mental  error,  i.e.,  it  is  Impossible  to  differentiate  between  the  two  models. 
Comparisons  of  the  overall  coefficient  model  with  the  pore  diffusion  control- 


5Marshall,  W.  R.,  Jr.,  and  R.  L.  Pigford,  "The  Application  of  Differential 
Equations  to  Chemical  Engineering  Problems",  University  of  Delaware 
Newark,  Delaware  (1947).  * 


ling  model  are  listed  in  Tables  3-5.  These  results  show  that  C/CQ  dif¬ 
fers  less  than  0.012  for  all  points,  and  for  most  points  the  difference  is 
much  less.  The  largest  difference  in  the  two  models  occurs  at  breakthrough; 
however,  this  difference  is  not  significant  since  the  experimentally  deter¬ 
mined  values  of  C/CQ  at  breakthrough  are  very  difficult  to  obtain  accurately 
due  to  the  very  low  concentrations  exiting  the  bed.  For  example,  with  a 
reference  peak  height  of  two  inches,  a  C/CQ  value  of  0.01  would  require  mea¬ 
suring  very  accurately  a  peak  height  of  0.02  inches.  The  range  of  UA  and 
shown  in  Tables  3-5  represent  extremes  not  found  in  this  study  but  are 
presented  here  to  illustrate  the  similarity  between  models.  The  comparison 
for  UA  =  5.9850  -  7.9785  with  =  15.53  is  representative  of  bolt  3  foam 
material,  whereas  the  comparison  of  models  for  UA  *  3.9900  -  5.9850  with 
Ka  =  7.76  -  11.64  best  represents  bolt  2  foam  material. 
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RESULTS 

Some  of  the  modeling  results  are  shown  in  Figures  2-6  for  bolt  3 
carbon  impregnated  foam  material.  The  model  prediction  is  shown  as  a 
solid  line  and  the  experimental  data  are  plotted  as  symbols.  Figures  2 
and  3  present  data  from  typical  one  and  three  layer  runs  on  a  conditioned 
sample  of  bolt  three  material  and  show  excellent  agreement  of  model  and 
experimental  results.  Figures  4,  5,  and  6  show  experimental  data  and  model 
prediction  for  three  types  of  conditions,  and  except  for  the  somewhat  scat¬ 
tered  data  in  Figure  6,  the  agreement  is  excellent.  These  results  indi¬ 
cate  that  the  model  is  adequate  for  predicting  breakthrough  curves  in  cases 
when  the  foam  material's  adsorptive  capacity  may  have  been  altered  by  sweat 
poisoning,  addition  of  water,  or  removal  of  water.  The  run  conditions  of 
Figure  6  represent  the  upper  range  of  flow  rate  and  concentration  used  in 
this  investigation,  and  the  experimental  scatter  is  probably  due  to  the 
very  rapid  rise  of  the  breakthrough  curve  with  limited  sampling  time.  Since 
and  the  rate  parameter  UA  are  calculated  from  the  experimental  data,  the 
poor  fit  by  the  model  in  Figure  6  must  be  due  to  inaccurate  data.  The  model 
is  helpful  even  with  scattered  data  as  it  fits  a  smooth  continuous  curve 
through  the  experimental  points. 

Figure  7  shows  the  effect  of  changing  the  equilibrium  adsorption  con¬ 
stant,  K^,  as  predicted  by  the  model.  The  characteristic  shape  of  the  curve 
is  the  same  but  the  area  is  different.  This  reflects  a  change  in  the  total 
amount  adsorbed  which  is  proportional  to  the  area  under  the  1  -  C/CQ  curve 


and  proportional  to 
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Figure  8  shows  the  effect  of  changing  the  overall  coefficient  by  +  25%. 
The  curve  for  UA  =  7.9785  goes  between  the  two  shown  and  is  left  out  to 
show  more  effectively  the  change  in  the  curve's  characteristic  shape.  The 
transmission  curve  is  not  very  sensitive  to  changes  in  UA  and  indicates 
that  small  variations  in  the  rate  parameter  due  to  experimental  error  may 
be  acceptable.  The  area  under  the  1  -  C/CQ  curve  is  the  same  for  both 
curves  drawn  and  reflects  the  fact  that  KA  (and  thus  total  amount  adsorbed) 
was  held  constant. 

The  three  independent  variables  that  influence  carbon  tetrachloride 
adsorption  are  temperature  (T),  concentration  of  the  carbon  tetrachloride 
in  the  inlet  gas  stream  (CQ)  and  flow  rate  through  the  sample  (Q).  A  cen¬ 
tral  composite  statistical  design  of  experiments  was  set  up  to  determine 
quantitatively  the  effects  of  these  variables  on  the  equilibrium  constant, 

Kft,  and  on  the  overall  mass  transfer  coefficient,  UA.  A  linear  fit  of  «A 
with  the  three  run  variables  showed  the  equilibrium  constant  to  be  a  signi¬ 
ficant  function  of  concentration  and  temperature  given  below  for  bolt  2 
carbon  impregnated  foam  material: 

Ka  =  17.4  -  0.58  CQ  -  0.12  T 

for  5.2  mg/1  £  C  <_  12.5  mg/1 

25°C  <  T  <  40° C 

with  a  standard  deviation  of  0.72  and  a  value  of  KA  (mean)  =  9.09. 

As  expected,  the  variation  with  flow  rate,  Q,  was  not  statistically 
significant.  These  values  of  KA  from  the  dynamic  adsorption  runs  (calculated 
from  the  zeroth  moment  equation)  compared  very  favorably  with  the  values  ob¬ 
tained  from  the  adsorption  isotherms  (from  McBain  balance).  However,  a 


slight  deviation  was  noticed  at  the  higher  concentrations  where  the  isotherm 
has  been  found  to  be  nonlinear.  This  is  due  to  the  assumption  of  a  linear 
rate  of  removal  made  in  the  adsorption  model,  that  is  not  valid  for  those 
concentrations  in  the  nonlinear  region  of  the  adsorption  isotherms. 

A  statistical  fit  for  the  overall  coefficient  for  the  statistical  de¬ 
sign  experiments  revealed  that  UA  was  a  function  of  1/CQ  and  Q: 

UA  =  2.10  +  4.80/Co+  0.75  Q 
for  5.2  mg/1  <_  Cq  <_  12.5  mg/1 

0.5  1/min  Q  <_  1.3  1/min 

with  a  standard  deviation  of  0.60  and  a  value  of  DA  (mean)  =  3.060.  The 
variation  of  UA  with  temperature,  T,  was  not  statistically  significant;  how¬ 
ever,  this  was  probably  due  to  the  small  temperature  range  (25°C  <_  T  _<  40°C) 
used. 

Figure  9  shows  the  results  of  the  above  correlation  produced  from  the 
statistical  analysis  of  38  runs,  and  the  data  points  shown  are  average  values 
of  UA  for  specific  run  conditions.  The  figure  shows  clearly  that  UA  is  in¬ 
versely  proportional  to  inlet  concentration  and  directly  proportional  to 
flow  through  the  sample.  This  further  indicates  that  the  overall  coeffi¬ 
cient  represents  more  than  one  resistance  controlling  since  DcA/R  should  not 
be  function  of  Q,  flow  through  the  sample,  and  k^A  should  not  be  a  function 
of  CQ,  inlet  concentration. 

The  statistical  analysis  for  bolt  2  foam  material  can  be  extended  to 
predict  the  parameters  UA  and  for  bolt  3  foam  material.  The  functional 
dependence  on  run  conditions  is  assumed  to  be  the  same  for  both  bolts  and 
the  intercept  used  in  the  empirical  correlations  is  determined  from  a  randomly 


picked  run  on  bolt  3  material.  Table  6  shows  the  result  of  breakthrough 
time  comparisons  for  bolt  3  material  based  on  the  predicted  values  of  the 
parameters,  and  except  for  one  run  at  a  low  concentration,  shows  an  excel¬ 
lent  agreement  of  predicted  and  measured  values  for  breakthrough  time. 


L  * 
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CONCLUSIONS 

While  the  use  of  an  overall  coefficient  model  is  not  limited  to  com¬ 
plex  adsorption  systems,  its  chief  advantage  is  the  simple  approach  to 
modeling  a  system  in  which  the  controlling  resistances  are  not  known  and 
are  not  easily  found.  The  definition  of  an  overall  coefficient  for  an 
adsorption  process  is  not  exact  in  its  description  of  the  process,  but  it 
does  provide  a  means  of  accounting  for  the  effects  of  the  three  individual 
resistances  in  consistent  units.  This  overall  coefficient  represents  ef¬ 
fectively  a  combined  resistance,  since  each  of  the  single  resistance  models 
fits  the  data  equally  well  when  the  rate  parameter  is  evaluated  from  the 
first  moment  equation. 

The  model  works  well  for  the  carbon  impregnated  foam  material  for  the 
range  of  flow,  temperature,  and  inlet  concentration  used  in  this  study. 
Extension  of  this  model  to  systems  using  an  adsorbent  bed  of  cloth  woven 
with  carbon  fibers,  beds  with  inerts  uniformly  mixed  in  with  an  unknown 
quantity  of  adsorbent,  and  similar  systems  should  be  no  problem.  This 
model  is  especially  advantageous  for  those  systems  in  which  flow  through 
the  bed  is  difficult  to  characterize,  the  amount  of  adsorbent  in  the  bed 
is  unknown,  and  in  which  various  physical  properties  (such  as  a,  e,  pp,  R) 
cannot  be  confidently  and  accurately  determined.  The  elimination  of  the 
need  for  moments  higher  than  the  first  moment  is  another  advantage  in 
treating  data  that  is  too  inaccurate  for  higher  (or  complete)  moment  analy¬ 
sis.  By  setting  up  a  simple  statistical  design  of  experiments,  one  can  de¬ 
termine  the  overall  coefficient,  UA,  as  a  function  of  the  various  run  con¬ 
ditions  and  with  this  information  the  dynamics  of  adsorption  can  be  charac¬ 
terized  mathematically. 
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List  of  Symbols 

A  surface  area  of  active  carbon  particles  per  unit  weight  of  foam 

2 

material ,  cm  /g 

2 

flow  area  of  foam  sample,  cm 

C  concentration  of  adsorbable  gas  in  the  interparticle  space,  mg/1 

C  .  concentration  of  adsorbed  gas  per  unit  weight  of  adsorbent, 

s  mg  CCl^/g  particle 

Ci  concentration  of  adsorbable  gas  in  the  intraparticle  space,  mg/1 
C  initial  concentration  of  CC1.,  mg/1 

0  H  o 

Dc  effective  interparticle  diffusion  coefficient,  cnr/sec 
kajs  adsorption  rate  constant,  g  particles/ (mgCCl^  -  min) 
adsorption  equilibrium  constant,  ml/g  foam  material 
Ka'  adsorption  equilibrium  constant  of  pure  carbon  particles,  ml/g 
mass  transfer  coefficient,  cm/sec 
L  number  of  layers  of  foam  material  in  the  bed 

M  zeroth  moment  of  t  vs.  {1  -  C/C  )  curve,  min. 

0  o  u 

M,  first  moment,  min^ 

3 

M0  second  moment,  min 

d  4 

M3  third  moment,  min 

Q  flow  through  sample,  1/min 

R  radius  of  spherical  particle  of  adsorbent,  cm 

r  length  coordinate  in  the  spherical  particle  of  adsorbent,  measured 
from  the  center  of  the  particle,  cm 

T  temperature  of  sample  and  gas,  °C 

t  time,  min. 

1/UA  overall  mass  transfer  coefficient,  ccCClA/(g  foam-min) 

V  interstitial  velocity,  cm/min. 

Wt  sample  weight,  g 

z  length  coordinate  . "  bed  of  adsorbent,  cm 
Greek  Letters 

a  cloth  porosity  (void  volume/total  volume) 

6  intraparticle  void  fraction  of  carbon  particle 

6  weight  of  particles  per  unit  weight  of  foam  material 


J 


void  fraction  representing  the  bed  volume  not  occupied  by  carbon 
particles  divided  by  the  total  volume 

apparent  density  of  foam  material,  g  solids/cc  total  volume 
apparent  particle  density,  g  carbon/cc  carbon 


Table  1.  Components  of  Experimental  Apparatus  (Figure  1) 

1.  Control  Valve  size  cy  =  0.038 

2.  Control  Valve  size  cy  =  0.00145 

3.  U-tube  Manometer 

4.  Carbon  Tetrachloride  Bubbler  Chilled  at  0°C 

5.  Orifice  Meter 

6.  Hook  Gage  Manometer 

7.  Control  Valve  size  cv  =  0.15 

8.  Rotameter 

9.  Rotameter 

10.  On-off  Valve 

11.  Rotameter 

12.  On-off  Valve 

13.  Thermocouple 

14.  Temperature  Equilibrating  Coil 

15.  Tangential  Entry  of  gas  into  Manifold  to  Facilitate  Mixing 

16.  Wire  Mesh  Obstruction  to  Improve  Mixing 

17.  Stainless  Steel  Manifold 

18.  On-off  Valve 

19.  Sample  Holder 

20.  U-tube  Manometer 

21.  Thermocouple 

22.  Rotameter 

23.  Control  Valve  Size  cy  =  0.15 

24.  Control  Valve  Size  cy  =  0.15 

25.  Plexiglas  Compartment 

26.  Exit  for  Reference  Sampling 


Table  2.  Comparison  of 


Table  3.  Comparison  of  Model  Results  -  Varying  UA 
Model  1  -  D  A/R  controlling  model 
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Table  6.  Model  Prediction  of  Breakthrough  Time  for  Bolt  3  Foam  Material. 
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FIGURES 


Figure  1.  Schema 


Figure  2.  Model  Prediction  and  Experimental  Data  for  On 


RUN  3-17 


Figure  3.  Model  Prediction  and  Experimental  Data  for  Three  Layer  Run. 


Figure  4.  Model  and  Experimental  Data  for  Different  Conditions. 


32. 5°C 


Figure  6.  Model  and  Experimental  Data  or  Different  Conditions. 


Figure  7.  Effect  of  Changing  K.  10  %  as  Predicted  by  Mathematical  Model 


Figure  8.  Effect  of  Changing  UA  25  %  as  Predicted  by  Mathematical  Model. 


Figure  9.  Correlation  of  UA  from  Statistica 
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Derivation  of  Zeroth  and  First  Moment  Equations 


mass  balance  of  the  adsorbable  component  in  the  gas  phase: 


...  3D  ,  3C . 

3C  3C  c  1— l  /  x . 

V  Tz  ~  Jt  W~  ol  3r  ‘ 


r=H 


(1) 


mass  balance  of  this  component  in  the  particle: 


D  dC .  3C . 

( - -  +  —  — L) 

6.2  r  3r 
3r 


dC  .  p  ac 

_ i_  ads  _  . 

at  ~  6  at 


rate  of  adsorption  (assumed  to  be  linear) 


(2) 


ac 


ads 


at 


k  (C.  -  C  .  /K  •) 
ads  l  ads  A 


(3) 


external  diffusion  boundary  condition: 
8C 

D  =  k  (C  -  C  ) 

c  dr  r=R  f  1 

internal  diffusion  boundary  condition: 


(**) 


ac . 

— r—  =0  at  r  =  0  for  t  >  0 
3r 


(5) 


and  the  initial  conditions 

C  =  0  at  z  >  0  for  t  =  0 


0^  =  0  at  r  >  0  for  t  =  0 


C  =  C  at  z  =  0  for  t  >  0 
o 


Derivation  of  Laplace  transforms  of  mass  balance  equations 
L[ C^( z ,r ,t) ]  s  C  j ( z ,r ,s) 


(6) 

(7) 

(8) 


D  d2C  dC.  P  __ 

Eq.  (2) — ^  -r  ( - -  *  -  —r—)  -  SC.  -  S  -f-  C  *0 

H  '  '  6,2  r  dr  l  6  ads 

dr 


(9) 


„  _  k  _  __ 

Eq.  (3)  “>  SC  .  =  k  ,  C.  -  -S2®.  C  ,  -A  C  , 

ads  ads  1  K, '  ads  '  ads 
A 


(..-sag—)  c 

S+k  ;  Ui 
ads 


substitute  (10)  into  (9) 


d  Cf  o  dc<  S0C.  Sp  k  , 

- L  +  —  — i. _ i.  _  — E.  (  „„  L  a£.®. . - )  c  =0 

p  r-  n  n  v  • '  <  u 


.2  r  dr  D  D  S+k  .  /K. 

dr  c  c  ads  A 

d2r 

bi  2  — 

=  - L  +  £  +  XC.  =0 

,2  r  l 

dr 


c  ads  A 


Transforming  boundary  conditions: 


3C.  dC 

(z,o,t)  =  0  (Z,o,S)  =  0 


ac.  k  _ 

•  ( z,R,S)  =  [cTs.z)  -  C.(S,R)] 
e 

Solution  of  (12)  is: 


(z,r,3)  =  E/r  sin  r  V~T 


-  ( z  ,R,S)  =  F,  [-^4  cos  R  VT  -  \  sin  R  VT] 
dr  R  r2 


Equating  (15)  and  (17)  “ ^ 


1c  , 

E  =  4-  C(C,,z)/-^4  nos  fc  VT  -  -4  sin  R  VT 


kf 

+  7r-?r  sin  R  VT 


(18) 
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Now  take  the  Lapalce  transform  of  Eq.  (l) 


,pr  3D  d C . 

_  V  —  -  SC*  -  — -£•  —  (— ±) 

dz  C  Ft  «  dr 


=  0 


(19) 


r=R 


from  (19)  and  (18) 


dC. 


dr 


1  (Z,R,S)  =  -  -f- - 


(k  /D  )(  VI/R  cos  RV^-  1/R2  sin  R\/T)c(S,z) 


\fT/R  cos  RV/^  **  1/R  sin  rVT  +  k  /D  R  sin  rV/T 

f  c 


(20) 


dC. 


dr 


“  ( z ,R,S)  =  MC(z,S) 


(21) 


where  M  =  1  - 


D  /k  VT  cot  r  VT  +  (1  -  D  /k  R) 
c  f  c  f 


(22) 


3k. 


Let  1/L  =  and  substitute  with  (21)  and  (22)  into  (19) 


.vk-sc-f  c  =  0 


§  *  <1  *  ?>  c  -  0 

Transform  boundary  conditions; 


(23) 

(2h) 


at  z  =  0,  c(0,t)  =  c 

solution  to  (2b)  is 
C, 


c(o,s)  = 


C(z,S)  -  ~  e  ~  S  2/V  -  M/L  2  =  ->e'S 


■M/Lz 


(25) 


(26) 


z/V  =  residence  time  in  bed,  assumed  to  be  very  small  when  compared 
with  M/Lz  so  that  (26)  becomes 
C 


C^S)  =  -f  e-M/L* 


(27) 
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c( z»t) 

Define  G(z,t)  =1 - - - 

o 

G(z,S)  =  i 

O 

Substituting  (27)  into  (29)  yields  the  transform  G(z  S) 


(28) 

(29) 


G(z,S)  =  |  (1  -  e"M/LZ: 


3D 


where  M/I,  =  ( i 


4  VT  cot  R  VT  +  (1  -  j^-) 


P  k 


1  -  -  r  *  s-rfffir:) 

c  ads  A1 


(30) 

(31) 

(32) 


Derivation  of  the  Oth  moment  of  G(S,z) 


using  Aris'  theorem: 

M  =  lim  G(z,S)  =  lim  l/S[l-e~M/Lz]  >  £ 
°  S-o  S-*-o  *  0 


applying  L'Hcpitals'  rule: 


M  =  lim  z 
0  S-hd 


d(M/L)  -M/Lz  _  _  d(M/L) 

— - s—  e  -  lim  z  •— — — 


dS 


S-*o 


dS 


expand  cot  P  V"T:  cot  x  =  . 

dlM_/Ll  _  3kf  1-c  dG(  A_)_  (  }  _ 

RV  (i  dS  h  r  '  1  “  D 


dS 


G(x)«  1 


t-2-  V”*  cot  R  \J1  +  1  - 


CfR 


^  *  1  -  X 


r  \n 


kfR 
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dGU) 

dS 


(1 


D  RA 
c 

3k  ' 


-2 

(-1) 


dA 

dS 


dA 

dS 


(-)  [f  ♦  (- 

c 


pk 


ads, 


k  .  /K.* 
ads  A 


(S  *  kadE/V)‘ 


M 

o 


_  lira  zd(M/L)  _  lira  l-£_  ,PcRA  .-2  V^_  r&_ 

~  S-K)  dS  '  S-K)  RV  a  1  k  3  '  3k,,. 

i  i  c 


P  k 

p  ads 


k«a.V 

<S*kMs*A? 


<6*VV>  (32> 

Derivation  of  1st  moment 

a,  dG(z  s)  /all  — M/liZ\ 

M  _  lim  *  _  lira  d(,l/S(l-» _ ) 

1  ”  “  S-K)  dS  "  "  S-K)  dS 


_  -  lim  d(M/L)  ze~M/LzS  -  (l)  (l-e~M/Lz)  v  0 
S-K)  dS  s2  0 


applying  L'Hopitals*  rule  twice: 


-  lim  d2(M/L)  ^o 
1  '  ^  '  dS2  '  2 


must  find 


-  lim  d2(M/L) 
3+0  dS2 


G(A) 


1  - 


D 


Z  SJO.  _  *%Y3L,  + !  _ 

r  3  ^ 
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0(a)  JS  1 


D  RA  R3D  A2 
c  c 


1- 


3k,  "  TjT 


,  .  .2 >  . .  ,,  2  2D  R3  2D  2R2 

lim  d  G(A)  _  lim  /dAv  r  c  c 

S^O  dg2  =  S-*0  "  W  LT5iT^  k  29 


)  - 


A  .  V 

ds2  kf3 


d2A  _  ,  s  Ppkads  ,  kads  ,  kads> ~3 

._2  "  K~d)  D  K.'  U  K  ' 

dS  c  A  A 


&  4  -  <-2»  V'a’V.* 

do 


lim  d2(M/L)  _  3kf  l-e_  lim  d2G(A) 


S"°  dS2 


RV  a  S-»0  2 

dS 


,,  ,,  2  0  +  p  K 

lim  /dA^  _  / _ p  A  , 

S-H)  'dS  '  D  ' 
c 


,  2 


substituting  the  above  three  equations  into  the  eq  lation  for  yields: 


2  2 

2  K  ^  ^  |y|  ~ 

«!  ■  *  »PV>?  BIT  *  (e  *  °pV>*  3^7  *  *  -I-  (33> 
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Explanation  of  Numerical  Work 

All  integrals  were  evaluated  numerically  on  a  computer  by  using 

Simpson's  rule.  The  integrals  listed  in  Table  (2)  with  an  infinite 

upper  limit  were  evaluated  until  the  exponential  damping  factor, 

exp(-a,),  reached  a  value  of  5  x  10  ^  and  a  subsequent  iteration 

-1* 

increased  the  area  by  less  than  5  x  10  percent.  Pigford  and  Marshall 
(6)  have  solved  the  set  of  equations  for  the  UA,  kadp ,  kf  model  and 
have  expressed  it  in  the  form  of  an  integral  with  a  finite  upper  limit. 
This  solution,  listed  below,  gave  the  same  values  for  C/CQ  as  the 
infinite  integral  expression  listed  in  Table  (2): 

f  e-"  I  (2(„,)1/2)d, 

O  J 

o 

hJKA'z 
V 


where  n 


.  .  .  1-t  .  UA 

•  1  =  t  h’  J  -  °p  ~  h  =  50 

r  A 


